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Recently, it was shown that in inclusive B — > decay, an angular decomposition provides 

three independent (g^ dependent) observables. A strategy was formulated to extract all measur- 
able Wilson coefficients in B — > Xa £^£~ from a few simple integrals of these observables in the 
low region. The experimental measurements in the low region require a cut on the hadronic 
invariant mass, which introduces a dependence on nonperturbative b quark distribution functions. 
The associated hadronic uncertainties could potentially limit the sensitivity of these decays to new 
physics. We compute the nonperturbative corrections to all three observables at leading and sub- 
leading order in the power expansion in Aqco/ nib- We find that the subleading power corrections 
give sizeable corrections, of order —5% to —10% depending on the observable and the precise value 
of the hadronic mass cut. They cause a shift of order -0.05 GeV^ to -0.1 GeV^ in the zero of the 
forward-backward asymmetry. 



I. INTRODUCTION 

The inclusive decay B — > Xs t^t' is highly sensitive 
to new physics, since it involves flavor-changing neutral- 
current interactions, which do not occur at tree level in 
the standard model (SM). It is described by the effective 
Hamiltonian 
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where Oi_6 are four-quark operators and 
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with Pl.r = (1 T 75)/2- Here we have neglected the 
s-quark mass. All short distance information, including 
possible new physics contributions, is encoded in the Wil- 
son coefficients, Ci. Thus, one can test the SM and search 
for new physics by extracting these Wilson coefficients. 
Two observables frequently studied for this task are the 
dilepton mass Up) spectrum and the forward-backward 
asymmetry [H, [2 [J ■ Recently, it was noted that a third 
observable, proportional to a different combination of 
Wilson coefficients, can be obtained from an angular de- 
composition and significantly increases the sensitivity to 
the different Wilson coefficients [3] . With the addition of 
this third observable, the precise measurement of the 
dependence becomes unnecessary. Instead, a few simple 
integrals of these observables are sufficient to deter- 
mine all measurable Wilson coefficients \n B ^ Xg £~^£~ 
with the data already available from the B factories. 

The presence of intermediate cc resonances, J/ip and 
^p' , restricts the portion of phase space that is amenable 



to a precise comparison between theory and experiment. 
There are two suitable regions, q^ < Wjy^ and q^ > m^i. 

The large q^ region is usually considered less favorable, 
because it has a smaller rate and suffers from large non- 
perturbative corrections. However, the experimental ef- 
ficiency is better there, and in Ref. [sj it was shown that 
by taking the ratio of the B ^ Xs £^£~ and B X^iP 
rates the nonperturbative corrections can be kept under 
control, so precise predictions are possible even at large 
q2. 

In this paper, we focus on the low q^ region, which ben- 
efits from a higher rate. The trade-off is that experimen- 
tally a cut on the hadronic invariant mass, mx < m'^*', is 
required to suppress the huge background from b c(— > 
s£^i')i^i>. This nix cut introduces hadronic uncertain- 
ties that can easily spoil the search for new physics in this 
decay. The problem is that the decay rate is put into a 
kinematic region where the usual local operator product 
expansion in powers of AqcD/'Tib is no longer applicable. 
Instead, the rate becomes sensitive to the motion of the 
b quark inside the B meson, which is described by non- 
perturbative b quark distribution functions (shape func- 
tions) [6,, 7}. (The large q^ region is unaffected by the 
mx cut.) 

The latest BABAR [i| and Belle P analyses use 
TO^J** = 1.8 GeV and m^J"* = 2.0 GeV, respectively. To 
obtain the total branching ratio in the low q^ region, the 
measurements are extrapolated to the full mx range us- 
ing signal Monte Carlo based on a Fermi motion model. 
The extrapolated measurements are routinely quoted to 
compare theory and experiment. This practice is ques- 
tionable because, as is well established in the context of 
inclusive B Xsj and B Xutv decays, any such 
extrapolation should not be considered reliable and can 
give at best a rough estimate of the effect of the mx cut. 

In the kinematic region of low q^ and small mx, one 
can calculate the inclusive decay rates using soft-collinear 
effective theory (SCET) [13, [ill- At leading order in 
AqcD/^ifc, they factorize into process-dependent hard 
functions ft,!*^! , a universal jet function J, and the uni- 
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versal soft shape function S |12l. Il3l|. i.e. 



(3) 



a result applied extensively in the study of inclusive 
B — > XulP and B — > decays. It was first applied 
to B ^ X, £+£- in Refs. HJ, [lH to study systematically 
the effect of the m'^^ on the spectrum and forward- 
backward asymmetry. In Ref. [15| it was shown that 
the cut on mx leads to a 10 — 30% reduction in the 
rate. This reduction is, to a good approximation, univer- 
sal among the different short distance contributions and 
one can take it into account accurately using experimen- 
tal information from B Xgj or B ^ Xuii', thereby 
maintaining the sensitivity to new physics. 

The largest irreducible hadronic uncertainties and 
universality breaking are expected to come from 
O ( Apon / mh) power corrections due to subleading shape 
functions |l6l.[r^.[T8'|. In this paper, we extend the anal- 
ysis of the three angular observables to incorporate non- 
perturbative shape-function effects arising from the mx 
cut, including the ©(AqcD/^^fc) subleading shape fimc- 
tions. 

In Sec. ini we briefly discuss the kinematics and the 
angular decomposition, defining the three observables 
HT,A,L{q^)- In Sec. Illli we discuss the separation of the 
perturbation series above and below the scale /i ~ mb, 
and our effective Wilson coefficients. In Sec. IIVI we 
present our results for Ht.a.l in the SCET region. The 
leading power contribution is given in Sec. IIV Al includ- 
ing the full NLL and partial NNLL perturbative correc- 
tions. The subleading power corrections are presented at 
tree level in Sec. IIVBI Their numerical impact is inves- 
tigated briefly in Sec. El and we conclude in Sec. IVII 



II. ANGULAR DECOMPOSITION AND 
KINEMATICS 
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FIG. 1: Phase space cuts relevant for B Xs£'^£~ in the 
Px plane. The measurements are performed in the orange 
(medium) region, where the mx and cuts overlap and 
Px <Px- 
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with q^ = q+q-. 

The functions Hi{q-^,p^) in Eq. ^ are independent of 
z, and are given by 
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The triple differential decay rate can be written as Q 



d^r 



dg2 dp+ dz 8 



= - [(1 + z^)HT{q^p^) + 2zHA{q\p-],) 



Here, 



+ 2(l-z2)i^4g^p+)]. (4) 

q^ = is the dilepton invariant mass. 

Ex T \px\, and z — cos 9. In or B" [S" or 
5+] decay, 6 is the angle between the ^+ [£^] and the B 
meson three-momenta in the £^ £~ center-of-mass frame. 
The q^ spectrum and forward-backward asymmetry are 
given by 
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The velocity of the B meson is — p^/mg- We 
define light-cone vectors n and n such that = = 
and p'^ = n ■ px, Px = n ■ px- For later convenience, we 
also define the leptonic light-cone variables 

q+^n-q^ rriB - Px , 



(8) 



In terms of the usual structure functions in the decom- 
position of the hadronic tensor. 
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= -g^^Wi + v^'v''W2 + ie^^a/sw^/Wa 

+ q^'q-'W^ + {v^q" + v"" q^')Wz , (9) 

the hadronic structure functions Wt,a,l in Eq. ([7]) are 
given by 



Wt = 4Wi, 

WA = -2{q+~q-)W3, 
WL^^q^Wi + [q+-q-fW2 



(10) 



Without any cuts, the phase space limits on q^ ^ p^, 
and z are 



< p^ < mB — yq^ ^ tiib 



'!< z < 1 



(11) 
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In the rest frame of the B meson, 



where 



2mB Ex = ml + m\ - , (12) 

so low corresponds to Ex ^ OiniB)- In conjmiction 
with the mx cut required by the experiments we have 
m\ <C E\ or equivalently pt^ <C ■ This is illustrated 
in Fig. [T] in the p^ plane. The measurements are done in 
the orange (medium) region, where the two cuts mx < 
2.0 GeV (dark red) and 1 GeV^ < < 6 GeV^ (light 
yellow) overlap. There, p^ ~ tob is large, while ^ 
1 GeV is small. This is precisely the kinematic region 
where shape function effects are important, as explained 
in the Introduction. 

More explicitly, a cut mx < m3^' corresponds to a 
dependent cut pj^ < p^^*" ^ where 
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The Hi with a cut on mx are thus given by 



H,iq';mT) = / dp+i^,(g^p+), 
Jo 

H,{qlql-mT) = /''dq^ HM^-mT) , (14) 

Jql 

where the phase space limit from Eq. (|lip is implicitly 
understood. 



III. SPLIT MATCHING AND EFFECTIVE 
WILSON COEFFICIENTS 

After the W, Z, and t are integrated out at a scale 
of order mw, the effective weak Hamiltonian in Eq. ([1]) 
is evolved down to the scale m^, where the decay rate 
is calculated by evaluating the matrix elements of the 
operators Oi. In this step, the contributions from the 
four-quark operators Oi_6 and Og can be absorbed into 
effective Wilson coefficients Cf^Q{q^) that are complex 
functions of q^ . In the evolution from mw down to 
7716, Cg{mi,) receives a ln(m^/m^) contribution from the 
mixing of O2, which formally enhances it to Cg{mb) ''^ 
0{l/as), while numerically |C9(TOb)| « Ciq. It is thus 
advantageous to separate the perturbation series above 
and below the scale mf,, such that below mi, the effective 
Wilson coefficients can be treated as 0(1) numbers. This 
is achieved by the "split matching" procedure introduced 
in Ref. in the context of matching on to SCET. 

The split matching can be thought of as first matching 
the effective weak Hamiltonian at a scale fiQ ^ mi, on to 
a sum of effective b ^ s £~^i~ currents, 
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In the second step, starting from Eq. the decay 

rate is calculated. In the local OPE treatment, the time- 
ordered products of the currents in Eq. (fTH)) are matched 
at the scale fii, ^ m^ on to a set of local operators, 
whereas in SCET, the currents are matched at fih ~ mt 
on to corresponding SCET currents, as we shall do in 
Sec. lIVl In either case, numerically one can take /Xf, = /io, 
while formally fxo and fib are independent scale parame- 
ters. For example, to estimate perturbative uncertainties 
they can and should be varied separately. 

While Jg is a conserved current and thus scale- 
invariant, the tensor current has an anomalous di- 
mension and is therefore taken to be at a fixed reference 
scale, fi = mf). To obtain a well-behaved perturbative se- 
ries, we use mf, in the IS" scheme [l^, although any other 
short distance 6-quark mass could be used instead. Since 
both sides of Eq. (|15p must be fio independent, and the 
currents are (by definition) fiQ independent, the match- 
ing coefficients (^^'"'^'(g^, /io) are also fiQ independent to 
the order in perturbation theory at which the match- 
ing is performed. Hence, the decay rate calculated from 
Eq. (jlSp is formally /ih independent and one can treat the 
^^-j 0(1) when counting powers of as below 
the scale pi,- This also means that we have to be careful 
with our terminology. As far as the Wilson coefficients 
are concerned, we stick to the usual _B — > £~^£~ count- 
ing, where, due to the formally leading l/a^ in Cgimi,), 
NNLL refers to 0{as). On the other hand, in SCET at 
fib and below, NNLL refers to the full two- loop 0{a^). 

Since the split matching happens at the level of cur- 
rents, it captures only finite virtual corrections, which are 
contained in the C™'^^{q^), and the universal IR diver- 
gent virtual and bremsstrahlung corrections, which are 
described by the currents Jyg. It does not incorporate 
finite bremsstrahlung corrections from operators other 
than 07,9,10, which must be added explicitly. In the local 
OPE their effect was shown to be small, around the 1% 
level [lOl^. In the SCET expansion they are both power 
and as suppressed and thus beyond the order we are 
working at. Similar considerations apply to electroweak 
corrections pTI. [23|. which are not included here. 

In Ref. [1| the C™'^'(g2,Mo) are decomposed as^ 



^ The full z dependence of these corrections, which may be known 
from the calculations of the authors of Refs. 20], has not been 
published, and so is not known for Hj'(q^) and Hj^{q^) sepa- 
rately, but only for HA{q'^) and the sum HtIq^) + H^iq^). 

2 In Ref. 3 the coefficients Cf^^^q^) are defined implicitly by 
absorbing all virtual contributions from Oi_6,8 into them and by 
requiring their /^o independence. That definition is equivalent to 
the one given here. The coefficients C™"'(g'^) in Refs. [l3 . ll5ll are 
equivalent to these except that C™'^((ji^) = {rag /mi,)C}p'^^{q'^). 
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Cr\q^^lo)=C, + F,{q^) + Goiq^), 
C[l^\q',^^o)^C^o, (17) 

such that all terms on the right-hand side of Eq. (|17ll are 
separately no independent to the order one is working at. 
The explicit expressions are collected in the Appendix of 
Ref. 4], and we do not repeat them here. To simplify 
our notation we suppress the /io dependence in the coef- 
ficients hereafter. 

The functions FT,g{q^) contain the virtual contribu- 
tions from Oi s and are known at NNLL order jH, 
[13, Hi,!!!, [Hill Hi] (up to small Oa.g contributions), 
while the GT^g{q^) contain nonperturbative Ojllmi) cor- 
rections involving the four-quark operators [30| ■ The lat- 
ter can be included in this simple form, but the final re- 
sults for the decay rates have to be re-expanded so that 
any terms of 0{as/m1, l/m^) are discarded. 

The coefficients Cy^g.io are real in the SM. They contain 
the dependence on the coefficients C7,9,io(At) in Eq. ([T]), 
i.e. 

C. = C.(Mo)^ + ---, 
nib 

Cg = Cg{flo) H , 

Cio = £"10 , (18) 

which are sensitive to new physics. On the other hand, 
the functions Fi{q^) and Gi{q^) are dominated by con- 
tributions from Oi^2 and thus are expected to be insen- 
sitive to new physics. Hence, the approach advocated in 
Ref. to search for new physics is to assume the SM 
everywhere and treat Cy^g.io as three unknown real pa- 
rameters to be extracted from data; it was shown that 
Hl(1,6), i/T(l,6), i7A(l,3.5), and 7/^(3.5,6) are suffi- 
cient for this purpose. This strategy has the advantage 
that the number of parameters is kept to a minimum and 
thus the sensitivity to new physics can be maximized. In 
addition, there is no dependence on a specific new physics 
model. New physics contributions will show up as incon- 
sistencies between the extracted values of Cr^g.io and their 
calculated SM values, or between overconstraining mea- 
surements (similar to the usual approach to overconstrain 
the Cabibbo-Kobayashi-Maskawa matrix) . 



IV. RESULTS 

In this section, we present our results for the three ob- 
servables iJr, Ha, and Hl defined in Eq. (g]) in the SCET 
region, <C pl^- We write their structure functions in 
Eq. dZl) as 

W,{q\px)^wl°\q\p+)^Wl^\q\p+) + --- , (19) 



where i = T,A,L and the superscript [n] denotes the 
order (AQco/^Tib)" in the power expansion. The leading- 
order wl^\q'^ , p~^) , involving the leading shape function, 
are discussed next, while the VF]^^(g^,pJ), containing the 
subleading shape functions, are discussed in Sec. lIVBl 



A. Leading order 

The leading-order structure functions factorize as 

wl'\q^pj,) = hf\q^p+,^^,) 

X Jdu;p^J{p^Lj,fj.i) S{pj^ - uj,fii) , (20) 

where p~ = rrib^q- = rrii, — q^ / {niB —Px) is the partonic 
light-cone momentum. The integration limits here and 
below are implicit in the support of the functions, which 
are nonzero only if their first argument is positive. We 
shall discuss each ingredient in Eq. ([^Dl) in turn. 

The hard functions hf\q'^,Px,^J.b) are different for 
each structure function. To obtain them, we start by 
matching the QCD currents in Eq. at the hard scale 
fif, ~ nib on to corresponding SCET currents, 

2=1,2,3 

j,^2rm ^,7(p-^^^)-^r.^^n^ (21) 

where Xn — W^^„ and — Y^hy are the standard 
collinear and heavy-quark fields in SCET, and p~ cor- 
responds to the large momentum label on the collinear 
quark field. We choose a slightly different set of minimal 
Dirac structures than usual, 

r^^, = Pr {iq,a''\q,{q-v^ - q^v")] . (22) 

The reason to use instead of n^^ for Egg is that it 
makes explicit the constraint from lepton current con- 
servation, which implies that for massless leptons only 
two coefficients contribute to the rate. For ^ there are 
only two independent coefficients from the start, because 
qtiJy =0 provides an additional constraint. 

The matching for general currents to 0{as) was car- 
ried out in Ref. [ll|. For the vector current, the two- 
loop 0{al) matching has become available only recently, 
through the work of several groups i32|, i33|, i3^] . We 
find 
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where the two-loop functions C^^\u) can be found in Ref. j34| . and as indicated they have to be evaluated at = 
For the tensor current, we find 

Clip ,Mb) = 1 Cf 
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The 0(0:3) corrections for the tensor current are not fully known at present, but since two-loop calculations for the 
vector current exist, the equivalent two-loop calculation for the tensor current should be feasible. From the two- loop 
computation of the ICVp terms in the & 37 rate [H, [3^, one can obtain the 0{a1) contribution to c[ at the 
point p~ = mb [stI. |38||. For the vector current, the a1 corrections to c^{p~,fib) for small or large p~ are to 
good approximation given by a constant shift. Assuming a similar behavior for the tensor current, we can obtain an 
approximate result for c[ at 0{al) in the low region. 



clip ,Mb)=Ci(p , ^Ib) + ^[hsimb, ^ib) - climb, ^J,bf] , 



(25) 



where cl is the result to ©(as) in Eq. ([24]), and hsimb,fib) is given to ©(a^) in Eq. (A4) of Ref [3§ |. 

The hard functions hf'^ are now computed by substituting the currents in Eq. (|2ip together with their prcfactors 
from Eq. (|15p into Eq. ([9]) for the hadronic tensor and factorizing out the matching coefficients times the trace of 

their Dirac structures. One then obtains wjul = h^l^l J S, with (writing for simplicity c}'^ = cf and rio,i = Fg^^) 



^[0] = ^ C-mcl* (jinc\ J2 cic] Tr 
a,b=7,9,10 i,i=l,2 



2 ^ a,i ^ ^ bj 



(26) 



The remaining matrix element gives the convolution of jet and shape function, J ® S. 
Taking the traces and the appropriate linear combinations from Eq. (|10p . we find [l^ 

h^°\q\p\,^,b) = CZ\q^) clip-, ^lb) + —Cr\q^) clip-, ^lb)\cl^ [c?(p-,Mb)]', 
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To evolve the coefficients from the hard scale fib ^ mb to 
the intermediate scale /i^ ~ mx ^ ^AqcDWb, we use 



hfiq\p+„^l,)^hl\q\p+^,flb)UHip-,^ib,^^^), (28) 



Ml. .-2 ^+ 



where the hard evolution factor [Tl[ sums logarithms be- 
tween the scales fib and fii and is known at NNLL. 

Next, we consider the convolution of jet and shape 
function in Eq. The jet function J(p~a;,/Zi) con- 

tains perturbative physics at the intermediate jet scale 
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fi^'^mx, and is known at 0{as) and ©(a^) [4l|. 

The leading shape function S{oj,n) is defined as'^ 



S{uj,fi) 



2mB 



{B\Oo{oJ,fi)\B) = {Oo{u:,fi))B, (29) 



where 



Oo(t^, A^) — by S{iD+ — S + uj)by 



(30) 



Here, by is the HQET b quark field, iD'^ is an ultrasoft 
covariant derivative, and 6 — ttib — rnb, so S(uj) has sup- 
port for Lj > 0. We use the full QCD B meson state \B) 
in Eq. which automatically absorbs into S{lu, /i) all 
subleading shape functions that would otherwise arise 
from time-ordered products of Oo(w,/^) with the power 
corrections in the HQET Lagrangian. 

The shape function contains both nonperturbative and 
perturbative physics. A method to combine all available 
perturbative and nonperturbative information was devel- 
oped recently in Ref. [11] . To do so, the shape function 
at the scale /i^ is written as 

S{tJ,fii)= duj' dkUs{uj - uj',fii,fj.A) 



X Oo(c^'-fc,MA)%), (31) 

where the hats indicate that the quantities are given in 
a renormalon-free short distance scheme. The function 
Co is perturbatively calculable at the soft scale ha, and 
is known at 0{as) MJMand 0{al) [13]. The soft 
evolution factor Us |40l 43l 44 1 sums logarithms between 
the soft scale fiA and the jet scale /i^. Finally, F{k) is a 
fully nonperturbative function, which can be constrained 
by data from B ^^7, B XylP and B XJP jH]. 

Combining the convolutions in Eqs. pT|) and ([^0]) . we 
define the perturbative function P {p~ , p'j^ , ^i) by 



P{P ,Px'/^») = J '^^ J duj' p J[p (pj - w),/ij] 

X Us{uj - uj' , fJ.i, ij-a) Co{uj' , ha) 
^5{k)+0{as), (32) 



and combining this with Eq. (|28|) we obtain 

X jdkP{p-,p+,-k,tJL,)F{k). (33) 



wf\q\p+,) = hf{q\p+,,^ib)UH{p-,^ib,^i^) 



With Eq. p7|) and the matching coefficients in Eqs. (1231) 

and (|25)) we have an approximate 0{a1) result for ft.^'^, 

which do not depend on While ft, I*' depends on C2, it 
has no soft photon pole and is thus completely dominated 



^ We use a different normalization of the \B) state from that in 
Ref. [H. 



by the vector current contributions, which are known at 
0{a1). An expHcit expression for P{p~ ,p~^, ^i) to 0{al) 
with NNLL summation, in any short distance scheme 
for the &-quark mass and kinetic-energy matrix element, 
has been derived in Ref. [1^. An explicit NNLL expres- 
sion for Uh can be found there as well. Hence, approxi- 
mate NNLL 0(0^) results are available at leading order 
in the SCET power expansion for all three observables 

Ht,a.l{i^,Px)- 



B. Subleading order 

At tree level and ©(AgcD/wb) six additional sublead- 
ing shape functions enter in the description oi B ^ Xu^v 
and B Xs-i [M EE 113, El, 111,113], and wiU 

also contribute to B ^ Xs £~^£~ . We refer to these as 
the primary subleading shape functions. The analog of 
the factorization theorem Eq. (1201) at ©(Aqco/'^h) was 
worked out explicitly in Ref. [47[. At ©(asAqcD/^^fc) 
an even larger number of additional shape functions ap- 
pears [13, E, [5l| . The split matching relies on the fact 
that for 07,9,10 we can treat q^ as 0(1) in the SCET 
expansion. If subleading contributions from other oper- 
ators are considered, it can be necessary to count q^ as 
parametrically small and to treat the photon as coUinear 
particle. In this case there will be additional four-quark 
operators with coUinear quarks coupling to the coUinear 
photon, giving rise to additional subleading four-quark 
shape functions [ij, [s^] • We shall restrict our discussion 
to tree level and the primary subleading shape functions. 

When we consider the ©(Aqcd/w;,) power corrections, 
the split matching is important for two reasons. First, it 
is convenient, because it allows us to think only about 
the two currents in Eq. ([TB)) . This implies that the fac- 
torization in Ref. [13] also applies to our case, and a large 
part of the results can be reused. More importantly, it 
provides us with a consistent way to work at tree level at 
the scale /ifc and below and neglect ©(asAgcD/wb) loop 
corrections in SCET, while at the same time keeping the 
full as corrections to the effective Wilson coefficients from 
scales fj,o and above, even when they multiply sublead- 
ing shape functions. In this way, we can avoid artificially 
large power corrections that arise simply from having to 
use different Wilson coefficients at ©(AqcD/wf,), and can 
instead use the same Wilson coefUcients at each order in 
the power counting. 

The calculation proceeds along the same lines as in 
the previous section, though here there are two sources 
of subleading corrections. First, the matching in Eq. (pij) 
now has to include subleading SCET currents. Secondly, 
when the time-ordered products are evaluated, there will 
be corrections from higher-order terms in the SCET La- 
grangian. Alternatively, working at tree level, we can 
directly match the time-ordered products of the effective 
currents on to the subleading shape function operators 
as in Ref. [13] ■ Of course, both approaches give the same 
results. 
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The operators that arise from subleading SCET cur- 
rents are 

O^'H ^^h{iD'',5{iD+-S + uj)}b,, 



Jduip J{p [p]^ - uj),^i){5 - uj)S{u) 

+ (35) 



P2{u) = A e^^, K &{:iD+ -S + oj)] 7^75 K . (34) ^^i^^ ^2(^) gives rise to a new subleading shape function, 



They come with the same jet function as the leading- 
order shape function. The contribution from 0^(a;) can 
be rewritten in terms of the leading-order result as 



dwp J(p {p\-u}),y.i)n^,{0'^{uj))E 



Aujp-J[p-{p+ ~ ij),fM) {P2{lo))b = F2{p+) + 0{as) . 

(36) 

The operators that are due to higher-order terms in 
the SCET Lagrangian are 



03{uji,LU2) = by S{iD+ ~S + uji) {iDi_f d{iD+ -S + LU2)bv, 
Pi'{uji,uj2) = ^bv S{iD+ ~6 + uji) ge±^xGf 6{iD+ - 5 + UJ2) 7t75 , 
0^^{uJi,uj2, Lo^) = [KT^ 5{iD+ -5 + LOi) 7^Pl5L] -5(19+ - 5 + L02) [K^i'^Pl S{iD+ -S + uj^) T^feJ , (37) 

where e^" = e^'^°'^nanf3/2 (with 60123 = 1), '^gG^^ = [iD^.W^], and s'l^ — (7^1^/4) Sus is an ultrasoft s quark field. 
These operators are associated with new jet functions that are known only at tree level. Combining their B matrix 
elements with their jet functions, we define 



F^{pjr) = J dui dtJ2 duj. 



Px - ^2 Px- ^1 
5{p+-uji) ^ 5{p+ - CJ2) 
Px Px - t^i 

47ra, r -r-r 1 



Fiip 



duji dw2 duj3 



27ri 

Anas 
27ri 



n 

= 1:2, 

n 

i=l:2,3 



(03(^1, Ws))^ , 
n^{P^{uJi,uj2))B , 

n 



■j = 1.2,3 Px - - j = 1^2,3 

1 



Px 



Px - - 



n 

J = l:2,3 



Px - ^3 + ie 



ni_,n,y{Oi^^ {uji,u!2,uj3))b , 



(.9M. + <.)(OrrK,^2,^3))B, (38) 



which correspond to the functions defined in Ref. [43|. There are also operators Pf , O2, P3, and O4, obtained from 
the above by interchanging the Dirac structure 1 7^75- They do not contribute because their matrix elements 
between B meson states vanish as a result of parity and/or time-reversal invariance. 

With the above definitions, we find the following ©(Aqcd/'ti;,) corrections to the structure functions: 



W^'\q',pi) 



wW(g2,p+) = -2CioRc|Cr'(9 
W^^\q^pj,)=[\Cr\q')\'+C!o 



[\Cr\q')\'+C!o 
4mb 



Fi{p+x) , Ft(p+) 



Re[C^=i((72)C^"=i*(q2)] 



mil p 
MPx) 



F,{p+) Ft{p+) 



^_x 

nib 



P' 



cl/„2^ 



F,{p+ ) , Ft{p+) 



+ ^Cr\q')^^^ 

q- p 



rrib p 
'qlF^ip+) qlFUpjc) + ^<l'F2ip+) 



rrib p 
Fi{p+x) , FL(pi) + 2P2(pi) 



mb 



P' 



(39) 



8 



Parameter 


Value 


ms 


5.279 GeV 


mt = [54] 


4.70 GeV 


Ai = Ai [38] 


-0.32 GeV^ 


A2 


0.12 GeV^ 


P2 M 


-0.065 GeV^ 


Ct 


-0.2611 


Cg 


4.207 


Cio 


-4.175 



TABLE I: Central values of input parameters. 

where we have used the abbreviations 

Fi{p+) = {p+ - 6)F{p+) + F2{p+) , 
MPx) 
Fdpjc) 



F3{p+)~F,{p+) + 2F^{p+) 
F3{p+)+F,{p+)^2F^{p+) 



(40) 



Note that in W^^ and w'^'^ only two different combina- 
tions of subleading shape functions appear, a property 
which can be exploited to construct particular combina- 
tions of observables in i? — > Xuiiy and B — > Xs£~^i~ for 
which the subleading shape functions drop out [s^. 



V. mx-CUT EFFECTS AT SUBLEADING 
ORDER 

In this section, we briefly investigate the numerical im- 
pact of the power corrections in Eq. ([39)l on the different 
observables, using the input values collected in Table[Il A 
detailed numerical analysis of the mx cut effects includ- 
ing an estimation of uncertainties is beyond the scope 
of this paper and is relegated to a dedicated publica- 
tion 

To obtain expressions for the leading and subleading 
shape functions, we follow the construction in Ref. [3^ . 
We give only a few relevant formulas here, and refer the 
reader to Ref. .38] for further details. The nonpertur- 
bative function F{k) entering the leading-order result 
Eq. ([55]) can be expanded as 



F{Xx) 



^ ^ fn (^) 



(41) 



where A ~ Aqcd is a free parameter and fn{x) form a 
complete set of orthonormal functions on [0,oo). We use 
the default value A = 0.8 GeV and /„(a;) from Eq. (48) 
of Ref. [38^. 

Since our main interest is in the corrections froni_sub- 
leading shape functions, we use a fixed model for F{k), 
obtained by truncating the series in Eq. ([1T|) at n < 2. 
For a given value of A, the remaining coefficients co,i,2 are 
determined by the 0th, 1st, and 2nd moments of F{k), 



I 



dkF{k) = 1, 



F(fc) 

--. F2{k)/mb 
F3{k)/mb 
Fi{k)/mb 




0.2 0.4 0.6 0.8 1 

k [GeV] 



1.2 1.4 1.6 



FIG. 2: Model functions used for the subleading shape func- 
tions F2{k) (dashed blue), Fi{k) (dotted green), and FA{k) 
(solid orange). The black solid line shows the model used for 
the leading-order function F{k). 



Ak k F{k) — rriB 



-.IS 



1S\2 



dfcfc"F(fc) = -^ + (ms-m^^*) 



(42) 



with ml given in the 15' scheme and X\ in the "invisible" 
scheme [38|. 

Very little is known about the subleading shape func- 
tions. Since the flavor of the light quark in the operator 
O^s does not match the flavor of the spectator quark in 
the B meson, we expect the functions F^ ^{p'^) to give 
only small corrections. Furthermore, since they arise only 
in combination with F3^4{p~\^) as in Eq. (|40p . we can as- 
sume that any small effect they may have will likely be 
washed out by the uncertainties in F3^4{p~^). We there- 
fore set FQfi{p~^) to zero in our numerical analysis. The 
first moments of the remaining functions are 

y"d/cf^2,3,4(fc) =0, 

J dkkF2,3Aik) - {-A2, 2Ai/3, A2} , 

'"dfc (fc - 5)2 f^2,3,4(fc) = {P2, 0, 0} . (43) 

For i^2,3,4(fc) we use a construction similar to Eq. (jiTI) . 

- 00 



F,{Xx) 



dx 



ri=0 



(44) 



which automatically incorporates the vanishing 0th mo- 
ment. The overall sign is determined by the sign of the 
first moment. To obtain a range of models for each func- 
tion we consider two cases, i and d\ 2 0, with all 
other coefficients set to zero. For each case, there are two 
solutions to the moment constraints from Eq. ()43p . pro- 
viding us with a total of four reasonably different models 
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1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 

m™* [GeV] m™* [GeV] 

FIG. 3: Relative corrections due to subleading shape functions as function of m^* for Ht{1, 6; m5^') (top left), -ffA(l, 3.5; m^*) 
and Ji"A(3.5,6;mx*) (top right), //l(1, 6; m^*) (bottom left), and r(l,6;m5'') (bottom right). The thin lines show the result 
of using the different subleading shape function models from Fig. [2] The thick black line in each case shows the center of the 
thin curves. 



for each function, which are shown in Fig. [2l When com- 
bined, these give 64 different sets of models for the sub- 
leading shape functions, which we use to illustrate their 
effects. We stress that the spread in the results obtained 
from these models should not be interpreted as a rigor- 
ous theoretical error, but merely as an indication of the 
rough size of the uncertainty expected from the unknown 
form of the subleading shape functions. A more detailed 
analysis will be presented in Ref. [ssj . 

To illustrate the effect of the power corrections, we con- 
sider their relative corrections to the lowest-order result. 



Hf\qlql-mT) 



(45) 



Here, Hi = {Ht,Ha,Hl,T} and the i/f^^' are ob- 



tained from Eqs. (|33p and (|39|) . respectively, correspond- 
ing to zeroth and first order in the power expansion. 
Since we consider at tree level only, we also use 

the tree-level result for i/]"' in the denominator for com- 
parison. Note that we keep the full NNLL expressions 
for the Wilson coefficients Cl'^^^iq'^) in both numera- 



tor and denominator (using the numerical expressions 
from Ref. [3). As already mentioned in Sec. IIV 51 this 
is consistent because of the split matching and is im- 
portant for maintaining the correct relative size of the 
different short distance contributions at ©(AqcD/'Ttb). 
Fig. [3] shows Ai7T(l,6;m™'), AFa(1, 3.5; m™*) and 
Ai7A(3.5,6;m5^'), Ai/L(l, 6; m5^*), and Ar(l, 6; m^J**). 
For Hl and F, the corrections are between 0% and 
— 10% with central values around —5% for to™* between 
1.8 GeV and 2.0 GeV. As expected from Ref. [H, the 
uncertainty in the correction increases for lower to™*. 
The corrections for Ht are somewhat larger with sim- 
ilar uncertainties. The reason is that in the combina- 
tion Frik) = Fz{k) — Fi{k) entering Ht the correc- 
tions from Fy,{k) and Fi{k) tend to add up, while in 
FL{k) — F^{k) -f F4(fc) entering Hl they tend to can- 
cel. 

Considering Ha, we see that the lower bin 77^1(1,3.5) 
receives a significant positive correction, above -1-10%, 
while the higher bin i?yi(3.5,6) receives only small neg- 
ative corrections of a few percent. The reason is that 
the C™'^^{q^)Cw term, whose absolute value decreases 
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3 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 



[GeV'] 

FIG. 4: The normalized function HA(q^;m''^^) for mT = 
2.0 GeV (solid and dark green) and mS"' = 1.8 GeV (dashed 
and light green) at lowest order (lines) and including power 
corrections (bands). 

(it becomes less negative), dominates Ha for very small 
q^. As increases, these corrections are compensated 
by a corresponding reduction of the C^^^{q'^)Cio term. 
This also results in a shift of the zero, qQ, where the 
forward-backward asymmetry dApB/dq"^ — {i/A)HA{q^) 
vanishes. In Fig. 2] we plot the ratio 

HA[q ,mx ) - , , (4b) 

r(l,6;m™^) 

i.e. HA{q^) normalized to the rate integrated over the low 
q^ region, as a function of for fixed *■ The black 
lines show the leading-order result ff]^'(<;^; m™') using 
m^J** = 2.0 GeV (solid) and m™' = 1.8 GeV (dashed). 
The green bands show the result obtained by including 
the subleading shape function corrections in both numer- 
ator and denominator, leading to a horizontal shift of 
about -0.05 GeV^ to -0.1 GeV^ with a similar uncer- 
tainty. This is the same size as the perturbative uncer- 
tainty usually quoted for q^. The size of the horizontal 
shift in the curve at q^ is not different from that at any 
other point in this q^ region. This is expected, because 
in the theoretical description of inclusive decays there 
is nothing special about the zero beyond the fact that 
HA{q^) happens to vanish there. 

VI. CONCLUSIONS 

In Ref. Q , it was demonstrated that the three observ- 
ables Hxiq'^), HA{q^), HL{q^) measured in the low q^ re- 



gion provide significantly better sensitivity to the differ- 
ent Wilson coefhcients than the rate dP/dg^ = HA{(f) + 
HL{q^) and forward-backward asymmetry dAps/dq"^ = 
{3/A)HA{q'^) alone. In the low q'^ region, the experimen- 
tally required cut on the hadronic invariant mass, mx, 
makes the measurements sensitive to nonperturbative b 
quark distribution functions, so-called shape functions. 
Rather than extrapolating the measurements to compare 
with theory, one should take the effect of the mx cut into 
account on the theory side. In this paper, we computed 
all three observables, Ht,a,l, in the low q^ region in the 
presence of an mx cut, including the leading and sub- 
leading shape function contributions. 

We used a split matching procedure to separate the 
perturbative corrections above and below the scale mf,. 
The perturbative corrections above mf, are taken into ac- 
count via uniquely defined effective Wilson coefficients 
C]"'='(g2) and are known at NNLL [0{as)] from the stan- 
dard calculation of i? — > Xs£'^i~ in the local OPE. Below 
the scale fib, the perturbative corrections at leading order 
in the power expansion are fully known at NLL [©(a^)] 
and approximately at NNLL [0(q;^)]. The subleading 
power corrections are included at tree level. 

While the effect of the mx cut at leading order can be 
taken into account model-independently by combining all 
constraints on the leading shape function from perturba- 
tion theory, together with available data from B Xgj 
and B Xu^c^v [s^, much less is known about the sub- 
leading shape functions, which represent a currently irre- 
ducible hadronic uncertainty. Depending on the observ- 
able and the value of the mx cut, the subleading shape 
functions induce corrections to the leading-order result 
of about —5% to —10% in the rates and a shift of about 
— O.OSGeV'^ to — 0.1 GeV^ in q^, with uncertainties of 
the same size. Hence, they must be accounted for to be 
able to obtain precise predictions for measurements of 
B ^ Xs i^l'^ in the low region. A detailed numerical 
analysis of the mx cut effects and their influence on the 
uncertainties in the extraction of the Wilson coefficients 
will be presented in a separate publication [55| . 
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